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The rich and complex optical response of high- 
temperature superconducting cuprates (HTS) is still at- 
tracting broad interest in the scientific community. In 
particular, much recent work, both experimental 1-2-3 ' 4 ' 5 
and theoretical^Ti?!?! 1 '. 1 ! 11 was focused on the temperature 
dependence of the Cu-0 plane infrared spectral weight 



W(Sl,T) = / dujatiuj). (1) 
Jo 

where <7\{u>) is the real part of the in-plane optical con- 
ductivity. In principle W(£l,T) depends on both the 
cut-off ft and the temperature T. When £1 — > oo, the 
f-sumrule for the optical conductivity is recovered. In 
this case W (x n, where n is the carrier density for a 
continuous system and therefore it is independent of T. 
In a tight-binding model W(Cl, T) is instead related to 
the carrier kinetic energyji 2 . and then retains a temper- 
ature dependence. In this framework, the behavior of 
W(Q,T) provides crucial information on the energetic 
balance at T c . Indeed, it has been reported that the 
onset of superconductivity both in underdoped and in 
optimally doped Bi 2 Sr 2 CaCu208+ a (BSCCO) is associ- 
ated with a gain in kinetic energyjiaSiSi as opposed to the 
slight kinetic energy loss in conventional BCS supercon- 
ductors. Moreover, restricted sum rules, where in Eq. 
l(T]l f2 < u>p, with io p the plasma edge, have been studied 
experimentally in BSCCO^ YBa 2 Cu 3 7 _ y (YBCO)^ 
and La 2 _ a; Sr ;I ,Cu04 (LSCO), 5 > both in the normal and in 
the superconducting phases. In the latter work, it has 
been shown that the quadratic T-dependence 

W(Cl, T) ~ W(n, 0) - B(fl)T 2 , (2) 

holds for any infrared cut-off SI between ~ O.leV and 
~ 1 eV, both in LSCO and in a conventional metal 
like gold. In both systems, when increases, the T- 
independent term obviously increases, while the thermal 
response B(fl) decreases. However, while in gold both 
W(Vl, 0) and B are controlled by the same energy scale 



(a hopping rate t), in LSCO the former scale is by one 
order of magnitude larger than the latter one. Moreover, 
in LSCO a sizable temperature dependence of W(Q,T) 
is found unexpectedly even for £1 ~ lj p , while in gold 
Bin ~ u p ) = 0. 

Expression J5J is easily recovered, at least for large cut- 
offs, if electron-electron interactions are neglected. If, for 
the sake of simplicity, one models the Cu-0 plane with an 
effective single-band tight-binding model with nearest- 
neighbor hopping t only, one has W(fl ~ u> p ,T) = 
— (Tre 2 )/(2d)Ekin, where e is the electric charge, d is 
the dimensionality of the system and Ekin is the car- 
rier kinetic energy. For a non-interacting system we can 
write Ekin as an integral over the density of states (DOS) 
and expand it a la Sommerfeld, obtaining Eq. J5J with 
W(uip) cx Ekin(T — 0). Under these simplifying assump- 
tions, Eki n (T = 0) oc t and B oc t^ 1 (for a quasi- flat 
DOS, e.g., B = 7r 2 e 2 /(24t d)), so that both quantities 
are controlled by the same energy scale t. As long as 
electronic correlations are neglected, this result cannot 
be significantly changed by more realistic models. It has 
been recently shown that the addition of a next-nearest 
neighbor hopping t' which controls the position of the 
Van Hove singularity cannot explain the observations in 
LSCO or BSCCOiS 

The above standard behavior of non-interacting sys- 
tems, in which the same t controls both W(ft,0) and B, 
is manifestly at variance with experimental evidence in 
cuprates. In particular, in LSCO^ Eki n {T — 0) leads 
to a to of the order of several hundreds meV, consistent 
with independent estimates of the bandwidth based on 
photoemission spectra or band-structure calculations. 13 
On the other hand, the experimentally determined value 
of B is unexpectedly large and leads to tx — 20 meVj 5 ^^ 
smaller than to by one order of magnitude. As it was 
proposed in Ref. |a, the existence of two distinct energy 
scales points towards the separation between low- and 
high-energy physics that is characteristic of strongly cor- 
related electronic systems, and suggests that the large 
difference between £o and tx can be due to the proximity 
to a Mott insulator. 
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The present paper is aimed at putting this hint on 
firmer grounds by calculating the optical conductivity of 
a strongly correlated single-band Hubbard model, by Dy- 
namical Mean-Field Theory (DMFT)ii We are aware 
that this approach will not provide a complete descrip- 
tion of the optical properties of cuprates, where spatial 
correlations and other effects definitely play a role. How- 
ever, we want to understand whether the Coulomb repul- 
sion, even in its simplest reliable description, provides the 
basic ingredient to explain the one-order-of-magnitude 
discrepancy between the T — spectral weight and the 
thermal coefficient B. We will assume reasonable values 
for the parameters of the calculation, but in this context 
it would be of little significance to finely tunc them in 
order to fit the experimental results. Those values, in- 
deed, should be certainly modified as other effects beyond 
DMFT were included. 

The Hamiltonian of the Hubbard model reads 



<ij><r i ' ' 

-//^(nii+njj), (3) 

i 

where Cj<r(ct ) are annihilation(creation) operators for 
fermions of spin a on site i, rii a = cl c^., t is the hop- 
ping amplitude, U is the local Hubbard repulsion, fi is 
the chemical potential and the sum in the first term is 
restricted to nearest-neighbors only. 

In order to solve iJSJ we resort to the DMFT, presently 
the most reliable tool to treat the physics of doped Mott 
insulators. DMFT is based on the neglect of spatial fluc- 
tuations, but fully retains local quantum dynamics, and 
becomes exact in the infinite coordination limit in which 
the self-energy is local, i.e., momentum- independent. 14 
Despite deliberately excluding phases with given spa- 
tial arrangement, DMFT allows for the simultaneous 
description of low- and high-energy features, and has 
provided the first unified treatment of the correlation- 
and density-driven Mott transitions. The DMFT maps 
the lattice model onto an impurity model subject to a 
self-consistency equation which contains the information 
about the original lattice. Here we consider, as it is of- 
ten done, an infinite-coordination Bethe lattice with a 
semi-circular DOS of half-bandwidth D. The impurity 
model has been solved by Exact Diagonalization at finite 
temperatures. The number of levels in the bath has been 
fixed to Nf, = 4 -=- 5, which, despite being small, gives a 
good description of thermodynamic properties, at least 
at intermediate temperatures^ Following Ref. ITEl we 
assume D = it = 1.2eV, and U = 3D = 12t, so that the 
antiferromagnetic superexchange J = At 2 /U ~ 100 meV 
is compatible with experimental estimates. As we discuss 
in the following, our choice is related to the message we 
aim to bring with this paper. 

Before turning to the analysis of optical properties, wc 
briefly remind how the strongly correlated phases aris- 



ing by doping a Mott insulator are characterized within 
DMFT. For large enough U/D, the half-filled system is a 
Mott insulator, with a spectrum made by two broad inco- 
herent features, the Hubbard bands, which are centered 
around ±Z7 /2 with a width of the order of the bare band- 
width. As the system is doped, a quasiparticle (QP) peak 
appears around the Fermi level, whose weight is measured 
by the residue Z = (1 — 3S(aj)/9w| aj=0 ) _1 7 be- 
ing the (momentum-independent) self-energy. Z, which 
tends to vanish when the Mott transition is approached, 
is small in the strongly correlated regime, reflecting the 
loss of metallic behavior due to the Coulomb repulsion. 
In the hole-doping case which we discuss here, the QP 
peak is located close to the lower Hubbard band. 

The knowledge of the single-particle spectral function 
can be directly used to compute the optical conductivity 
(and the associated spectral weight) because vertex cor- 
rections vanish in the infinite coordination limit in which 
DMFT becomes exact^ and ai(u>) (summed over all the 
spatial directions) is simply given by 



cr^uj) = — J de J dv N(e)V(e)Im[G(e,v)] 

x MG{e, V + u)] f{v) - f{V + W) (4) 

UJ 

where N(e) is the DOS of the non- interacting system, and 
V(e) = (D 2 — e 2 )/3 is the square current vertex of the 
infinite-coordination Bethe lattice we are considering, 16 
and G(e, oj) = (uj — e + fi — E(w)) -1 is the retarded Green 
functions of the lattice system. 

The behavior of crj.(a>), already discussed in Refs. 
IrlTl clearly mirrors the above described spectral fea- 
tures of a doped Mott insulator, and presents three main 
features, namely (i) a Drude contribution at low frequen- 
cies (say up to 0.1.D) arising from optical excitations 
within the QP resonance, (ii) a Mid Infrared (MIR) peak 
(at uj ~ 0.5 -r 1.0-D) resulting from the transitions be- 
tween the lower Hubbard band and the unoccupied part 
of the quasi-particle peak and, finally, (iii) an insulating- 
like contribution at higher frequencies (at u) ~ U) mainly 
due to transitions between the Hubbard bands. 

Now we discuss how the above described behavior re- 
flects in the temperature dependence of the optical spec- 
tral weight W(fi, T). As it has been done for experimen- 
tal data, we consider two values of the cut-off: f2 = 0.1D, 
which is expected to include in the integral only the con- 
tribution from the Drude peak, and £1 = 1.5D, which is in 
a sense the plasma edge of our model, i.e., contains both 
the Drude weight and the MIR bands, excluding from the 
integral the small but finite high-energy contributions. 

In Fig. [T]the ratio W(T)/W(Q) is plotted for both the 
above f2 values and for different doping levels. One can 
see that our model calculation nicely follows the T 2 be- 
havior in the range 0.008D < T < 0.022Z? (i.e., between 
100 and 300K), as experimentally found in cuprates»i*2i&& 
The agreement with observations is not limited to the 
T 2 law, but extends to two important properties which 
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FIG. 1: W(Q., T) as a function of T 2 for the Hubbard model, 
for Q. — 0.1D (left panel) and £1 = 1.5D (right panel), nor- 
malized to its zero temperature extrapolation. Symbols are 
the results of the DMFT calculations, lines are best fits to 
them. 
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TABLE I: Results of the fits in Fig. 1 for selected doping 
values. First and second columns: values of the coefficient B 
(in units of 1/-D) for low- and high-frequency cut-off, respec- 
tively (note that, for U = 0, B ~ 2/D in the whole doping 
range here considered). Third and fourth columns: corre- 
sponding spectral weights. Last two columns: kinetic energy 
and quasiparticle weight. 



are not accounted for by non- interacting calculations. 
Namely, we find that (i) the inclusion of strong cor- 
relations determines a substantial increase in the coef- 
ficient B controlling the thermal excitations, and that 
(ii) the enhancement is more pronounced for the lower 
value of Q. The results for B(Q = 0.1D, 1.5D) and 
W (n) = W(Q,T = 0) determined by fitting the nu- 
merical data with Eq. are reported in Table I. For 
the sake of comparison we also report the T = values of 
the quasiparticle weight Z and for the full kinetic energy 
E kin ,le., W{Q^oo,T = Q). 

The large value of B reproduced by the DMFT can be 
mainly ascribed to the renormalized effective bandwidth 
2D* = 2ZD determined by correlations. 23 If however 
the same effective bandwidth controlled also the zero- 
temperature spectral weight in our calculation, a strong 
inconsistency between theory and experiment would ap- 
pear for this last quantity, whose experimental value is 
compatible with photoemission estimates of the hopping 
parameter, in turn of the order of the bare hopping. 

The explicit calculation of Wo(f2 = 1.5D) is instead 
in very good agreement with the experiments, as shown 
in Fig. [21 Therein the experimental points are obtained 



FIG. 2: Zero-temperature spectral weight for Q — 1.5D as a 
function of doping, for both the strongly correlated Hubbard 
model (dots) and for the observations in LSCO (circles). The 
latter data are obtained by integrating the a\{u) of LSCO in 
Ref. E3 

by integrating the (J\(uj) of LSCO reported for different 
doping in Ref. 0. Fig. shows that, even if the coher- 
ent hopping scale is significantly reduced by correlation, 
Wo is much less affected and it is still large due to the 
contribution of the MIR band. Indeed it involves inco- 
herent hopping processes which contribute to the average 
kinetic energy, even if they have no role in QP transport. 

It is worth to notice that theoretical and experimental 
points clearly follow a similar behavior as a function of 
doping, with a decreasing weight when the Mott insulator 
is approached. 

In this regard, we expect that the inclusion of spa- 
tial correlations beyond DMFT and of the antiferromag- 
netic superexchange J would result in a non vanishing 
renormalized bandwidth, as predicted, e.g., by slave bo- 
son mean- field theories i^a At the same time, once the 
DMFT approximation is relaxed, vertex corrections have 
to be considered in the optical conductivity. The fact 
that in Fig. [^J the calculated values are systematically 
lower than the experimental ones reflects our choice to 
fix U / D without adjusting its value to data. Anyway, 
the ability of the present DMFT calculation to give the 
correct order of magnitude of Wq suggest that the effect 
of the superexchange and the inclusion of vertex correc- 
tion may partially cancel each other, as it indeed happens 
within zero-temperature slave boson approaches^ 

In the final part of this paper we like to consider a 
quantity suitable to measure the variation with temper- 
ature of the optical spectral weight at large f2. Denot- 
ing here as W{T) the spectral weight W(fl,T) where 
£1 = Lo p for the experiments and f2 ~ 1.5D for the 
calculation, we plot in Fig. the relative variation of 
W(T) between T = and T = 300 K AW/W = 
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FIG. 3: The relative variation of the spectral weight be- 
tween the lowest T and 300 K, as a function of doping 
for various cuprates (full symbols) is compared with present 
DMFT calculations (squares) and with the predictions of non- 
interacting models (tight binding model (boh line), flat band 
(boh-boh line)). The dashed line is a guide to the eye. The 
simple inclusion of correlation effects allows one to reproduce 
the observed absolute values with no need f or fitt ing parame- 
ters. Data for LSCO are obtained from Refs. H ill , for BSCCO 
from Ref. (BSCCO i) and from Ref. Q (BSSCO2). 

[W(300K) - W(0)]/W(0). In the same Figure, we report 
for comparison the prediction of a non interac ting near- 
est neighbor model (dotted line, see also Ref. 10). Fig. 
13 clearly shows that: i) such ratio has similar values in 
different families of cuprates, indicating that this quan- 
tity is essentially material-independent; ii) the thermal 
change of W does not depend on details of the Fermi 
surface, which are instead dependent on both material 
and doping; iii) DMFT, with the introduction of strong 
correlation effects, reconciles the calculations with the 
experimental infrared behavior of the cuprates, which in 



no way can be reproduced by more conventional hopping 
models. 

In conclusion, we have shown that the inclusion of 
strong correlation effects, already in the simplest the- 
oretical scheme given by the Hubbard model, removes 
the inconsistency of non-interacting models and natu- 
rally accounts for the bifurcation between a T = spec- 
tral weight which becomes rapidly comparable with the 
non-interacting value, and the coefficient of its tempera- 
ture dependence where a renormalized effective hopping 
plays the major role. This neither means that the Hub- 
bard model in the DMFT approximation can account for 
all properties of cuprates, nor that our analysis has the 
ambition to be quantitative. However we believe that a 
non-perturbative treatment of correlations, such as the 
DMFT approach, is a basic building block for a the- 
ory of the optical response and of the transport prop- 
erties of high-T c materials. In our view, the next step 
would be the inclusion of spatial correlations by means 
of non-local extensions of the DMFT such as Dynamical 
Cluster Approximation (DCA)2i or Cellular DMFTM In 
the latter approach, the effective bandwidth is no longer 
bound to vanish at the Mott transition, as it should 
also depend on the antiferromagnetic exchange coupling 
J. As mentioned above, we finally expect that the ef- 
fect of the band broadening induced by J will be par- 
tially compensated by vertex corrections to the optical 
conductivity&ii These improvements would ultimately 
lead to a better understanding of the role played by 
charge and spin modulations in the optical response of 
cuprates. 
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